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Abstract — We consider the broadcast relay channel (BRC), 
where a single source transmits to multiple destinations with 
the help of a relay, in the limit of a large bandwidth. We address 
the problem of optimal relay positioning and power allocations 
at source and relay, to maximize the multicast rate from source 
to all destinations. To solve such a network planning problem, 
we develop a three-faceted approach based on an underlying 
information theoretic model, computational geometric aspects, 
and network optimization tools. Firstly, assuming superposition 
coding and frequency division between the source and the 
relay, the information theoretic framework yields a hypergraph 
model of the wideband BRC, which captures the dependency 
of achievable rate-tuples on the network topology. As the relay 
position varies, so does the set of hyperarcs constituting the 
hypergraph, rendering the combinatorial nature of optimization 
problem. We show that the convex hull C of all nodes in the 
2-D plane can be divided into disjoint regions corresponding 
to distinct hyperarcs sets. These sets are obtained by superim- 
posing all k-th order Voronoi tessellation of C. We propose an 
easy and efficient algorithm to compute all hyperarc sets, and 
prove they are polynomially bounded. Then, we circumvent the 
combinatorial nature of the problem by introducing continuous 
switch functions, that allows adapting the network hypergraph 
in a continuous manner. Using this switched hypergraph ap- 
proach, we model the original problem as a continuous yet non- 
convex network optimization program. Ultimately, availing on 
the techniques of geometric programming and p-norm surrogate 
approximation, we derive a good convex approximation. We 
provide a detailed characterization of the problem for collinearly 
located destinations, and then give a generalization for arbitrarily 
located destinations. Finally, we show strong gains for the optimal 
relay positioning compared to seemingly interesting positions. 

Index Terms — Low SNR, computational geometry, network 
optimization. 

I. INTRODUCTION 
Next-generation wireless standards, such as 3GPP Long 

Term Evolution- Advanced (LTE-A) standard HI, propose re- 
lays as a mean to extend cellular coverage or to increase data 
rates. More specifically, LTE-A defines relays of Type I as 
coverage-extension relays which allow a base station (BTS) 
to reach uncovered users in a cell, and relays of Type II as 
relays which allow to increase the communication rate of a 
user already covered through a direct link to the BTS 0), 
0. In terms of cellular deployment, a natural and practical 
question arises as to where the relay node should be deployed. 

ITMANET - 6915101: This material is based upon work under subcontract: 
18870740-37362-C, issued by Stanford University and supported by the 
Defense Advanced Research Projects Agency (DARPA). 



In this paper, with the downlink of a cellular system with 
relays in mind, we address the aforementioned question for 
the broadcast relay channel (BRC), which consists of a single 
source broadcasting to multiple destinations with the help 
of a relay. In this paper, we focus on the wideband regime 
of wireless relay networks, also denominated low signal-to- 
noise ratio (SNR) regime because power is shared among a 
large number of degrees of freedom, making the average SNR 
per degree of freedom low. We would like to point out that 
addressing the low-SNR regime is relevant in next generation 
cellular systems. Indeed, considering LTE, large bandwidths — 
up to 20 MHz — can be supported by all terminals ( |3]). Due 
to power constraints in the low SNR regime, relays appear as 
a meaningful and natural way to increase rate and reliability. 

Previous results on wireless systems in the low-SNR regime 

include the capacity of point-to-point additive white Gaussian 

noise (AWGN) channel J4|, and multipath fading channel J5]- 

ifTTl . both equal to the received SNR: C Fading = Cawgn = 

p 

-£- = lim VFlogfl H ); the capacity of the mul- 

N <> W^oo V WNq ' 

tiple input multiple output (MIMO) channel [12), 03); the 
capacity region of the AWGN broadcast channel (BC) |fl4|- 
02), and AWGN multiple access channel (MAC) 07); and 
bounds on the capacity of the non-coherent multipath fading 
relay channel 08). From these works, a conclusion can be 
drawn on wireless systems in the low-SNR regime; the major 
impairment in the low-SNR regime is neither multipath fading 
nor interference, but noise, which is in contrast with the high- 
SNR regime. Formulating the argument more concretely, in 
the presence of multipath fading in the low-SNR, the same 
rates as the AWGN system with the same received SNR can 
be achieved using non-coherent peaky signals whereas spread- 
spectrum signals perform poorly. Moreover, the low-SNR 
regime is not interference-limited: in particular, all sources in 
the low-SNR MAC can achieve their interference-free point- 
to-point capacity to the destination. Based on this observation, 
the authors proposed in a recent work |fl9) an equivalent 
hypergraph model for the low-SNR AWGN MAC and BC. 
Then they used these models to build an achievable hyper- 
graph model for a more complex wireless network with fixed 
sources, relays, and destinations, and showed that optimizing 
power for maximizing multiple session rates boils down to a 
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straightforward linear program. 

In this paper, we take a step forward by simultaneously 
optimizing the relay location and the power allocation, to 
maximize the multicast rate from a source s to a set of 
destinations T. Using concepts from information theory, com- 
putational geometry and network optimization, we develop a 
comprehensive and efficient way to solve this problem, which 
can be broadly divided into three parts: 

1) BRC hypergraph model: We propose a hypergraph 
model for the low-SNR BRC, which depends on the 
topology of the network, essentially the placement of 
nodes on a 2-D plane. Given the source and destina- 
tions positions, computing the hyperarcs in the BRC 
hypergraph model requires to get the ordering of nodes 
in increasing distances from the source and relay, for 
all relay positions (as the relay is not initially given). 
This problem can be modeled as an ordered k-nearest 
neighbor problem, for which we propose a solution based 
on superimposing the Voronoi tessellations of all k — 1 
orders, where k spans the destination set. 

2) Continuous hypergraph variations: For fixed source 
and destination positions, when the relay position varies, 
the the network hypergraph changes accordingly ren- 
dering the problem combinatorial. Consequently, tradi- 
tional network optimization techniques cannot be applied 
directly, as they assume a fixed given hypergraph. To 
circumvent this hard combinatorial nature, we introduce 
continuous switch functions which allow to change the 
network hypergraph in a continuous manner as the relay 
position changes. Ultimately, this allows us to cast the 
problem as a continuous optimization problem. 

3) Convex approximation: The resulting continuous net- 
work optimization problem is non-convex. However, 
using geometric programming (GP) and p-norm approx- 
imation techniques, we provide a good convex approxi- 
mation of the original problem to which standard convex 
optimization techniques can be applied lEUl . It should be 
noted that the problem is NP-Hard in its original form 
mainly due to combinatorial nature and continuous non- 
convex constraints. 

Hereafter, the paper is as follows. In Sections HI] and 
III, we build the system model and formulate the general 
problem, respectively. In Section [TV] we solve the problem for 
collinearly located destination nodes, and introduce algorithms 
to compute distinct hypergraphs for various relay positions. 
Section [V] extends to the general problem case for an arbitrary 
topology, finally leading to the conclusion in Section [VI] 

II. LOW SNR SYSTEM MODEL 

Notations: IN and R denote the sets of non-negative integers, 
and real numbers, respectively. Let m £ IN and M m = 
{1, . . . , m} and Let S be a set, the indicator function of S 
is defined by l s (x) = 1 if x € S, t s {x) = if x £ S. 

In this section, we first recall the equivalent hypergraph 
models of the wideband BC and MAC, then we use them to 
build an achievable hypergraph model of the BRC, and finally 



we formulate the optimization problem. In hypergraph models, 
a hyperarc [u, V\v^ ■ ■ ■ Vk) of size k connects a transmitter u 
to an ordered (increasing order of distance from transmitter 
u) set of k receivers {vi, Ua, ■ ■ • , Vk}, all of which can decode 
a message sent over the hyperarc equally reliably. Here, u ^ 
{vi, i»2, • ■ • , Vk}- Two hyperarcs are disjoint if either they have 
different sources, or different ordered receiver sets or both, e.g. 
(u,ViV2) is disjoint from (u, V2V1) and {u' ,v\V2). Messages 
sent over any pair of disjoint hyperarcs are independent. A 
hyperarc is said to be activated if its capacity is non-zero. 

A. Wideband BC and MAC model 

Earlier results on multiple user channels show that BC and 
MAC are not impaired by interference in the low SNR regime. 

1) Equivalent hypergraph of the wideband AWGN BC, 
1 14 /- ^76l/ . K2l\l : Superposition coding is known to achieve 
the capacity region of the AWGN BC. In the wideband limit, 
the rates achieved by superposition coding boils down to the 
time-sharing rates, rendering time-sharing as optimal. 

Consider the BC channel with source s, two destinations 
T = {d\, CZ2 } in Figure [1(a)] D uv > is the distance of node 
v from u, and let a <E [0, 1] be the power-sharing factor at 
source. Then both destinations can receive the common rate 
amm{hi, h^}^, and the most reliable destination can also 
receive a bonus private rate (1 — a) max{^, where h 

is the path loss factor, A*o the channel noise and P is the total 
source power. This motivates the equivalent hypergraph model 
of the wideband AWGN BC in Figure [1(b)] The wideband 
BC hypergraph model contains three hyperarcs: the common 
hyperarc from the s to d\ and ^2 with capacity equal to the 
common rate R c , a private edge from s to d-y with capacity 
equal to the private rate R\ (if d\ is more reliable than c?2, 
i.e. h\ > h\, and to otherwise), and finally a private edge 
s to c?2 with a capacity equal to the bonus rate R2 if c?2 is 
more reliable than d\. Note that the two private hyperarcs 
(associated with rates R\ and R2) cannot exist simultaneously 
(as either D s \ < D S 2 or D S 2 < D s i): thanks to the indicator 
functions in the capacity expressions, only one of the private 
hyperarcs can be activated for a given topology. In the general 
case of a wideband AWGN BC with n destinations: 

• For an arbitrary unknown topology, the full hypergraph 
model contains 2™ — 1 hyperarcs, from the source to every 
possible subset of destinations. 

• For a given known topology, only a subset of these 
hyperarcs are activated simultaneously. Indeed, a given 
topology yields a given ordering of the destination set in 
increasing order of reliability. Consequently, only n hy- 
peracs are simultaneously activated for a given topology: 
one private arc of size 1 to the most reliable destination 
i.e. (s,d\), and one common hyperarc of size k to the 
k most reliable destinations for all k £ {2, ...,ra} i.e. 
(s,di..d k ). 

2) Equivalent hypergraph of the wideband AWGN MAC 
KTTj : Consider two sources S\ and S2 and a single destination 
d in Figure |l(c)| In the wideband regime, the large number 
of degrees of freedom renders negligible interference, and 



(a) BC (b) Wideband BC (c) MAC (d) Wideband MAC (e) BRC (f) Wideband BRC 

equivalent hypergraph equivalent hypergraph achievable hypergraph 

Fig. 1. Wideband Multiple User Channels. The BC rates are i?i = (1 - a)h\ j^t] h 2 t+00 [{h\ ), R2 = (1 - a) -^1 [0h 2 y{h\), R c = a min{h\ , h\}j^. 
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allows all sources to achieve their point-to-point capacity to 
the destination, as with frequency division multiple access 
(FDMA). Thus, the respective capacities of s x and s 2 are 
Ci = h\j± and C 2 = h\j^. As shown in Figure [T(d)l the 
equivalent hypergraph model contains only two edges, one 
from s\ to d with capacity C\ and one from s 2 to d with 
capacity C 2 . In the general wideband MAC with n sources, 
the hypergraph model consists of n hyperarcs of size 1 with 
non-zero capacity, from each source to the destination. 

B. Wideband BRC model 

Consider the broadcast relay channel in Figure 1 (e), where 
a source s transmits to a set of n destinations T = {dijigw 
with the help of a relay r. We assume that all nodes are 
equipped with a single antenna. The source and the relay 
have given respective average power constraints P s and P r , 
and they transmit in two different frequency bands, W s and 
W r respectively, so as to respect the half-duplex constraints at 
the relay. During each time slot, s transmits a new codeword 
which is received by the relay r relay and T (all destinations); 
r processes the signal received from s in the previous time 
slot and retransmits it to T; the destinations use the signals 
they received directly from s and through r to decode a new 
codeword. 

The wireless link between two nodes u 6 {s,r} and v £ 
{r, di, . . . , d n } is modeled by an AWGN channel. In other 
words, when node u transmits a signal x u (t), node v receives 
a signal y vu (t) = h vu x u (t) + z v (t) where h vu = is 
an attenuation coefficient modeling pathloss, and z v {t) is a 
white Gaussian noise process with power spectral density N$. 
Note that although we consider AWGN channels, the low- 
SNR analysis could be extended to multipath fading channels: 
indeed, it was shown in lfl8ll that in the wideband multipath 
fading relay channel, the same rates can be achieved as in 
the wideband AWGN relay channel with the same average 
received SNR on each link. 

The AWGN BRC consists of two BC components in series 
for s and r, the BC from s to {r, d\ , . . . , d n } in red in 
Figure |l(e)| and the BC from r to {d\, . . . , d n } in blue — 
and of n parallel MAC components, such as the MAC from 
{s,r} to g?i represented by the sum of the red and blue lines 
resulting into a green line. 

As in 1191 . we make the assumption that the source s and the 
relay r are constrained to transmit using the scheme that would 



be optimal for their respective wideband BC-components: 
s transmits using superposition coding in the source band 
W s ; r decodes the messages it received from s, and then 
retransmits using superposition coding in the relay band W r \ 
each destination di decodes by using the interference-free 
signals it received from s and r. 

Under these constraints on the communication scheme, 
the resulting hypergraph model of the BRC fl9l is simply 
the concatenation of the equivalent hypergraphs of the BC- 
components and the MAC components. We will denote this 
hypergraph as Q(Af, *H), where, JV D T = {di, .., d n } and H 
is the set of hyperarcs. The set T~L can be partitioned into two 
disjoint sets: one is H s of source hyperarcs emanating from 
s, and the set H r of relay hyperarcs emanating from the r, 
where T~L S U W r = %. Figure |l(f)| illustrates the hypergraph 
in the case of a given topology with two destinations. In this 
figure, we assume that r is the closest node to s, followed by 
di and then d 2 , and we show only the activated hyperarcs. 

It should be pointed out that this BRC hypergraph model 
is only an achievable model, and not an equivalent model. 
Indeed, in the case of a single destination, the BRC boils down 
to the relay channel, and it was shown in ifTSl that with a 
different coding scheme, it is possible to achieve a higher rate 
in the wideband relay channel than any rate obtained by the 
scheme in fl9l . Thus the BRC hypergraph model proposed in 
this paper provides only an achievable rate region, but not the 
full rate region of the BRC. However the relaying scheme in 
1 19], and the associated hypergraph model, have the benefit to 
easily extend to large complex network. 

III. General problem structure 

Given a topology of the set of nodes Af\r, and the 
aforementioned achievable hypergraph model of the wideband 
BRC, we recall the following question: What is the optimal 
relay position and power allocations at s and r, that maximize 
the multicast rate R m from s to the destination set T? Here, 
the multicast rate is the rate experienced by the least reliable 
destination in the set T, and is given by its min-cut. To 
solve the problem in full generality, we propose a two-stage 
approach, as follows: 

• Pre-processing: The pre-processing stage computes all 
distinct H s and H r , respectively, for all positions of the 
relay inside the region being considered on the 2-D plane, 
given by hyperarc sets H s and H r , respectively. Since, 
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only a subset of these hyperarcs are active when the 
relay is in a certain region, we associate each hyperarc 
(u, V) G H, along with a continuous switch function 
f u v- The switch function f u y activates the hyperarc 
(u, V) by taking the value 1 when it should exist and 
deactivates the hyperarc by taking value when it should 
not exist. In this section and section IV, we devise 
efficient algorithms to compute all the distinct hyperarcs. 
Once the hyperarc set H is constructed, we can then 
compute all the possible paths from the source to each 
destination. The total number of paths from s to a 
destination di € T will be denoted Ki, and the rates 
on these paths will be written {r*, ...,r^.}. 
• Optimization: The second stage involves solving a net- 
work flow optimization problem. After obtaining %, the 
multicast rate maximization problem for G(N, H) can be 
formulated as: 

Program (A): maximize (R m ) 
subject to: R m < r i7 Mi G 1N„, (1) 

n<J2 r k' VielN n , (2) 
fc=i 

max r* < y uV , i G lf„, V(u, V) G U, (3) 

((i,ki)\ki£(u,V), 1 

fcie[i,*r<]) 

Vuv < c u vUv, V(u, V) e H, (4) 
where, c uV G C uV , M(u, V) G H. (5) 

Here ([TJi implies R m is the minimum among the total 
rates experienced at all destinations, (ffji says that the rate 
for destination di is the sum of rates on all the paths from 
s to di. (0) captures the network coding constraint, and 
the switch function f u y in activates and deactivates 
hyperarc (u, V) as the optimization algorithm goes from 
one relay position to the other to maximize R m . (O 
implies that the capacity of hyperarc (u, V) G H is 
determined by implicit constraints of power and distance. 
Note that, it is because of continuous switch functions that we 
have a continuous optimization problem, else we would need 
combinatorial constraints to capture the right hyperarcs that 
need to be activated for each relay position. 

In the sequel of this section, we describe the two stages in 
detail. At first, we show that the optimal relay position lies in 
the convex hull C of nodes JV\r = {s, T}. 

A. Convex hull lemma 

We prove hereunder an intuitive lemma which helps in 
formulating the problem in a geometric sense. 

Lemma 1: Given Q(J\f,T-L), the relay location that maxi- 
mizes the multicast rate R m from s to T lies inside the convex 
hull of the nodes Af\r. 

Proof: We prove this lemma by building an intuitive 
argument. To start with, consider the four node system of 
{s, r, d%, d%} in Figure [2] The only node that is allowed to take 
its desired location is the relay r. The convex hull C of the node 
set Af\r = {s, di, c^} is shown as the shaded area in Figure |2] 
Consider the arbitrary position outside C in the above scene 




Fig. 2. Three line line segments joining (s, d\,d,2) form the convex hull 
C (blue region). The hyperarcs (s,r), (r,di) and (r, d\d,2) are shown with 
dashed arrows. 

e.g. r is placed outside C. The hyperarc sets for this position 
of r are given by H s = {(s, r), (s, r, di), (s, r, d\, d 2 )} and 
H r = {(r, rfi), (r, di, cfe)}, where T-L, assuming that the 
ordered sets of nodes in increasing distances from s and r 
are given by (s, r, dx,d2) and (r,d\,d2), respectively. 

Consider the rates on the path {(s, r), (r, di)} from s to di, 
which are given by: 

where, P s and P r are source and relay powers. It is clear 
that by moving r towards the boundary of C, i.e. line segment 
joining s and di, total rate on this path (min(i? sr , R r i)) could 
be increased. The triangle's inequality corroborates this fact, 

D sr + D rdl >D sdl , (7) 

(Rsr + Rrl)\{D sr +D rdl >D 3dl } 

< (R sr + Rrl)\{D sr +D rdl =D 3dl }- 

© is the triangles inequality for the triangle Asrd! and this 
implies (O, which states that the rate from s to d\ on the path 
{(s,r), (r, di)} could be increased by simply bringing the r to 
towards the line segment joining s and d\ . It is straightforward 
to see that this also increases the rate for all other receivers 
in system (consider triangles As,,^). 

This reasoning can easily be generalized to any arbitrary 
position of the relay outside the convex hull C, and to any ar- 
bitrary number of destinations \T\ > 2. Thus, we conclude that 
for any given instance of BRC, the relay location maximizing 
the multicast rate lies inside or on the border of the convex 
hull of {s + T} nodes, but never outside. Hence, proved. ■ 

Lemma [TJ implies that we only need to consider relay 
locations in the convex hull C of nodes Af\r, to maximize the 
multicast rate. There are efficient algorithms for constructing 
a convex hull of n + 1 points, c.f. ll22ll and references within. 

B. Pre-processing algorithms 

The pre-processing stage consists of three sub-stages. For 
all relay positions in C, the first sub-stage computes the source 
hyperarc set TL S , the second sub-stage computes the relay 
hyperarc set Hr, and the last sub-stage computes all the 
source-destination paths Vd; G T. Now, we develop algorithms 
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Fig. 3. (a): T = {di , efo , } ordered set T with three concentric circles 
ci for each destination, (b): Shows the ordered set of nodes w.r.t. r for the 
disc and two rings. Starting from the disc, the H s is computed and then for 
each ring the two new hyperarcs are added. 

to compute T~L S and H r and give upper-bounds on the number 
of distinct hyperarcs in H s and W r . 

1) Source hyperarcs (H s ): We first prove a lemma on \H S \. 

Lemma 2: The number of distinct source hyperarcs inside 
the convex hull C is upper bounded by 3n— 1, where n = \T\. 

Proof: Consider a BRC with s and an ordered set of 
three destinations T — {di,d^,ds} (\T\ = n = 3). Let 
be the circle centered at s passing though di. An example is 
illustrated in Figure [3ja), the three circles partition the 2-D 
plane into rings and discs that are given by disc C\ and two 
concentric rings 72-21 and 72.32. 

For the positions of r inside these areas there are dis- 
tinct sets of source hyperarcs. Computing the source hy- 
perarcs when r is inside these regions and starting with 
the disc c\, we simply get a set of 4 hyperarcs H s = 
{(s, r), (s, rdi), (s, rd\d 2 ), (s, rdi cfeda)}. Each time r crosses 
the border of circle Cj and enters the ring TZu-i, there 
are two hyperarcs that change and hence the new hyper- 
arcs must be added, i.e. {(s, d±), (s, dir)} for 72.21 and 
{(s, d±d 2 ), (s, did 2 r)} for 72.32, (refer Figure |3b)). This is 
due to the fact that, when r enters a new region, the ordered 
set of nodes in increasing distances from s of the new region 
is different from the previous region in only two places, e.g. 
for the disc C\ and ring 72-21 the ordered sets are given by 
{s 1 r 1 d\ 1 d 2 ,dz) and (s,di,r,d 2 ,d 3 ), respectively. 

Thus, the maximum number of distinct source hyperarcs 
that can exist for all relay positions in C is given by (n + 1) + 
2(n— 1) = 3rt — 1. Hence, proved. ■ 

At this point we would like to highlight a couple subtleties: 

(a) If some d t € T are equidistant from s, then their respective 
circles coincide, hence reducing the number of disjoint 
rings, and the number of distinct hyperarcs becomes less 
than 3n — 1. Thus, Lemma [2] is an upper bound. 

(b) The simple algorithm outlined in Lemma [2] builds all the 
possible hyperarcs efficiently in the sense that only distinct 
hyperarcs are added to H s along with their respective 
switch functions, thus avoiding any redundancy. 

Switch Function: The activation/deactivation of a hyperarc 
can be performed by the switch function. For instance, the 
switch function associated with hyperarc (s,rdi) (Figure [3} 
is fsri = (1 + 7e ( ~ 72srl) )~ 1 , where z sr i = D s i - D sr and 
7 >> 1. When r is in C%, z sr \ is positive, and f sr \ ~ 1, 




Fig. 4. The 4 destinations are shown with the cells showing ordered pairs (in 
increasing distance) in the plane after superimposing 1st, 2nd and 3rd order 
Voronoi diagrams. The dashed polygon forms C. 

thus hyperarc (s, rd\) is active. Similarly, when r is outside 
Ci, z sr \ is negative, and f sr i ~ 0, thus hyperarc (s,rdi) is 
deactivated. Notice, that for the hyperarcs inside a ring (e.g. 
TZu-i), the switch function will be a product of two functions, 
each for the regions of concentric circles that make this ring. 

2) Relay hyperarcs (T-L r ): The relay hyperarcs are deter- 
mined by the ordering of node set T with increasing distance 
from r. Thus, we need to partition the 2-D plane into disjoint 
regions where the ordering of T with respect to r stays the 
same. This is equivalent to computing the order-Zc Voronoi 
tessellations of C for all k 6 !N n _i, and then superimposing 
them to get the ordered set T of nearest neighbors in C, l23l . 
It is known that the superimposition of Voronoi tessellations 
results in convex disjoint areas (polygons in our case). Once 
the ordered set T for each disjoint region is computed, the n 
activated relay hyperarcs for each region are obtained the same 
way as in Section lTl-Al for the case of a BC with a given known 
topology. Thus, the algorithm outlined in Lemma 1 could be 
used to generate H r - Figure 2J illustrates the superimposed 
disjoint regions of ordered destinations with respect to r for 
n = 4. The simplest way to compute these regions (ref. 
[24 1) is to draw the perpendicular bisector of every destination 
pair (di,dj) in T 2 . This method has the complexity of order 
0(n(n—k)), where k £ IN n -i- Hence, we obtain the partitions 
of C for distinct ordering of the set T with respect to r from 
which we can generate H r - 

The pre-processing in almost all network planning problems 
is computationally heavy as there is plenty of time up-front 
compared to real-time applications. In our case it fits better as 
all the computations are of polynomial order. 

In the next sections, we formulate the problem of optimal 
relay positioning as a non-convex network flow optimization 
problem and provide a good convex approximation. For sim- 
plicity and clarity in understanding, we divide the problem 
into two cases. The first case is for collinearly located des- 
tination nodes and the second case is for arbitrarily located 
destinations. 
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IV. COLLINEAR CASE 

In this section we develop the method to solve the sim- 
pler version of the problem where the destination nodes are 
collinearly located. It helps understanding the main concepts 
and underlying algorithms, and ultimately leads to the solution 
for the arbitrary case. 

A. Pre-processing 

1) Convex hull C: For the collinearly located destination 
set T, the set T could be ordered in the increasing order of 
abscissa, for instance with the left most node being at the 
origin. Hereafter, we assume that the left-most (respectively 
rightmost) node di G T (d r G T) is situated at the origin 
(respectively right most at horizontal axis), and the rest of the 
destination nodes are on the positive x-axis, and finally that 
the source is in the positive quadrant. Note that, s could be 
the leftmost node compared to any di G T, in this case s 
could be assumed to be on the positive vertical axis and the 
set T would accordingly be placed on positive horizontal axis 
with di not being at origin. Since all di G T are collinearly 
located, the convex hull will always be the triangle A s d lt i r (ref. 
Figure |5(a)) . Thus, C is always given by only three inequalities 
in this case. 

2) Source hyperarcs: As explained in Section Ull-Bll the 
source hyperarcs are functions of the source-destination dis- 
tances. Consequently, the algorithm outlined in Lemma [2] 
could be used to compute H. s . 

3) Relay hyperarcs: Also, as shown in Section IHI-B2I we 
compute the perpendicular bisectors of every destination pair 
(di,dj) G T 2 , to compute the superimposed convex disjoint 
regions of ordered destination sets w.r.t. r (ref. Figure |5(b)) . 
For the collinearly located nodes, the computation of the set 
T-L r is greatly simplified due to parallelism of all bisectors. 
The following Lemma is just an easy and straightforward 
formalization. 

Lemma 3: For \T\ = n collinear destinations, the total 
number of distinct relay hyperarcs in C is upper bounded by 
n 2 = n + 2(j3 - 1), where j3 = (™) + 1 is the number of 
bisected regions. 

Proof: Given a set T of n colinear destinations, the 
maximum number of bisectors are given by (")■ Then the 
total number of bisected regions is given by (3 = (™) + 1, 
as shown in Figure |5(b)| Since, crossing each bisector only 
changes two nodes in the ordered destination set, using the 
algorithm outlined in Lemma\2\(ref. Figure|3]l, we can compute 
all distinct relay hyperarcs to T, \H r \ — n + 2(/3 — 1) = n 2 . 
Hence, proved. ■ 

With only little more formalization of Lemma [3] we can 
device easy algorithms for computing the set W r . Note, that 
the switch function for each relay hyperarc can be computed 
in a similar manner as for the source hyperarcs. At the same 
time, a particular switch function could constitute two sub- 
switch functions each for two perpendicular bisectors. 

4) Source-destination paths: After successfully computing 
all distinct source (H s ) and relay (TL r ) hyperarcs for a given 
system, we now need to compute all the paths from ,s to all 




(a) Source hypergraph. 




(b) Relay hyperarcs. 




(c) Source-destination paths. 

Fig. 5. Pre-processing. The triangle A s( j 2C j shows C, with circles and 
perpendicular bisectors dividing C is closed and disjoint sets, (a): Shows 
ordered 4-tuple set for each region carved by the circles c\ , C2 and C3 as 
Ci, 7^21 1 and 7£4oo, respectively, (b): Shows the ordered 3-tuple sets 
of destination nodes with respect to r. (c): Shows the previous 2 figures 
superimposed showing the disjoint convex regions with the ordered sets 
of closest nodes with respect to s and r, respectively. The 2-tuple (C, r) 
represents the ordered sets for each region, respectively. Here, n,r2,r , 3 and 
T4 represents the four regions in (b). 
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destinations di G T, in order to successfully cast our problem 

as a network optimization program. We prefer a path based 

formulation as opposed to a more basic and standard link based 

formulation because the path based formulation is far well 

suited for convex approximations of originally non-convex 

network optimization programs, in our framework. 

There are many efficient ways (polynomial time algorithms) 

to compute the paths from the set H s and H r . For simplicity, 

we prefer here to take all combinations of the hyperarcs in 

H s and H r . Not all these paths will be active for a certain 

relay position, but the switch functions will take care of 

activating/deactivating the paths. In this way, we get an upper 

bound (3n — l)(n 2 ) = \H S \ X \H r \ on the paths from s to 

each di G T, out of which only a certain number of paths 

will have non-zero min-cut, i.e. activated hyperarcs. We define, 

fl = [1, (3n — l)(n 2 )] as the set of all paths from s to every 

di G T. This makes the problem size bigger, but saves cost 

of activated path computation. Note that the total number of 

paths are polynomially bounded in our model. 

B. Optimization. 

In this section we formalize the problem of optimal relay 

location maximizing the multicast rate from s to T. The 
optimization constraints can be grouped into two categories: 
the posynomial constraints that can be easily rewritten using 
exponential transformation as convex constraints (Geometric 
Programming); and the non-posynomial constraints that can 
only be approximated as convex constraints. Since, all the 
constraints are coupled through variables, almost all the vari- 
ables and hence constraints will go through the exponential 
transformation. Now, we classify and discuss the troubling 
constraints in our formulation. 

1) Hyperarc rate constraints (Posynomials): We show an 
example of hyperarc rate constraints. Consider the hyperarc 
(s, d\rd2) in the scenario in Figure [5(a)] which is active when 
the relay is inside the ring 72.21- The non-convex rate inequality 
can be expressed as: 

Vs,\r2 < -f^T^rfslr2 (9) 



D%N ' 

where f slr2 = / s \ r2 / s 2 lr2 , / s \ r2 < (l + 

, 7e (-7^ 2 )r 1 ) Dsr _ Dsl = g i 



fslr2 < 1 



and 



z sir2 — Ds2 — D sr - Notice, how fli r2 and / 2 lr2 and their 
respective z variables are different. When, / s 1 lr2 and / 2 lr2 
take the value as 1, the hyperarc will have a non-zero min- 
cut. Rewriting them together, 



ft lr 2 /sir 2 



< 1, 



Rslr2D? 2 No 



< 1, 



(10) 



Pslr2f}i r 2fslr2 

/iir2(l+7e ( " 7 * ilrt) ) <1, (ID 
/ s V 2 (l+7e("^) 



D sr D s i < z s i r2 , 



< ft. 



< 1, 



(12) 
(13) 



Note, that inequalities ([Tol l. (fTTT i. ( fT2l and (fT3l are posyno- 
mials, and D s i is a constant. Using GP transformation these 
inequalities can be easily converted to convex constraints. Sim- 
ilar argument goes for switch functions of all other hyperarcs. 



2) Distance function constraints (Non-posynomials): Vari- 
ables D uv in rate inequalities represents distance functions, 
given by e.g. 



s) 2 + (Vr ~ Vsf 



D M 



(14) 



where, (x s ,y s ) are fixed coordinates of s and (x r ,y r ) are 
the variable coordinates of r. The negative coefficients in ( fl4b 
prohibits the use of GP techniques. There are techniques to get 
around this problem ll25l . Il26l . but the extra pre-processing 
cost incurred is very high in addition to the introduction of 
many new variables and combinatorial constraints. 

We prefer to handle the issue in a simpler manner by 
approximation. Let, the only variable transformed using GP 
in ( fT4b be D sr . Then, we can rewrite 



u lr + v sr ^ e 2D "~, x r - x s < u sr , y r - y s < v s 



(15) 



The first inequality in dT3T > is non-convex. Using the p-norm 
surrogate approximation ( fl27l ) for (fTBT l, we get 



u 2 +v 2 NP 

*~ar ~ sr 
2D' 



< 1, 



(16) 



where p £ [1, +oo). Over a compact set of variables and in the 
limit of p — > oo, (TToT l becomes convex. In our case, for values 
of p = 4 or 5, we get good approximation. Note, that only the 
first inequality in (fT~5b needs to be approximated, and since the 
variables (u sr , v sr , x r , y r ) don't undergo GP transformation, 
the rest of the inequalities in ( fT5l ) remain linear. 

All other constraints in the program are posynomials, as we 
will shortly see, so they can easily be transformed into convex 
constraints. It should be noted, that it is only because of the 
use of switch functions that the program becomes continuous. 
In addition, carefully designing the switch function results in 
posynomial hyperarc rate constraints. 

3) Network Optimization problem formulation: Since there 
are ft = (3n — l)(n 2 ) number of paths for each destination 
di G T only a subset of them will actually be active (i.e. 
with min-cut > 0). as the rate on path q to destination di, 
where q G O. Recall, T~L = T-L s U % r and the total rate to a 
destination di be defined as r, = J2 q en r g i - Also, let v m G V 
be the farthest node from i for hyperarc (u, V) G H. Then, 
the optimization program is, 



Program (B): maximize (i? m ) 
subject to: R m < rj,Vi G W n , 



(17) 
(18) 



max Rl < 



((<,?)!(«, V)e 5 ) 9 ' D$ V _N 



f uV ,V(u,V)£H, (19) 



fuv<flvfuv, V(u,V)eH, (20) 
fLv < (l + ie-b**^)- 1 , l€ [1,2], V(u, V) G H,(2l) 
z l uV < D UVl - D ur , I G [1,2], (u, V) G H, (22) 
u 2 uV + v 2 uV < D 2 ur , V(u,V)eH, (23) 
Puv <P S , p w<Pr, (24) 



((u,v)eHs) 



«U,V)£H r ) 



g 



where, x r — x u < u u y, V(u, V) g H, (25) 
Vr - Vu < v u v, V(m, V) g H, (26) 
i r > 0, y r > 0, y r < Xx r , y r + X'x r < r\. (27) 

In the above program (x r ,y r ) are variable relay coordinates 
and (27) captures constraints that make C. Program (B) is a 
non-convex program expressed in posynomial and polynomial 
inequalities. Applying GP transformation to the following vari- 
ables {r^, P u y, D u y, fuv}> p-norm approximation to con- 
straints ( 1231 and leaving the rest of the variables unchanged, 
we get the following convex approximation, 

Program (C): maximize min r'^ 

d ' eT \qen 

Noe (r': v +<xD' UVm -P^ v -f> uV ) < 1 

subject to: ' (28) 

\/q g (u,V),V(u,V) g U, 

e ULv-fuv-fuv) < i ) V(u,V) g (29) 
e (/iV) + e (/iV-74v) < i ? / e [l,2],V(ti,V) g H, (30) 

V(!i,y)6H, (31) 



4j + 

2 

uV 



"7 



J uV 



E 

((«,v)eH„) 



e p ^ < P„ 



< 1, V(u,V) <E H, (32) 



((«,V)e-H s ) 



where, we have used the GP transformation x' = log (a;) (a; is 
the original variable of program (B)) on certain variables. 

Program (C) is a convex approximation of program (B) with 
no underlying combinatorial hard structure. The approxima- 
tion is only coming from constraint type d32b using p-norm 
surrogation technique that gives a convex approximation to the 
constraint d23l in program (B). 

Note, that the objective function is modified in (C), instead 

of having the sum of positive exponential terms, we have 

replaced it by a sum of linear functions, which is far easier to 

maximize. The maximizers of the function maxtxi +x?) also 

x 

maximizes the function max(e Xl + e X2 ), over a compact set 
X with certain particular characteristics. This is generally not 
true, but in our case due to network coding constraints with 
certain simple tricks it can be proven that it holds true. Due 
to the lack of space we omit the detailed proof. 

As we know, with the increasing value of p, the program 
approaches a complete convexity with zero duality gap, thus 
standard convex optimization algorithms could be used to 
solve problem (C) with increasing accuracy. The optimal 
values of program (B) could be easily constructed from the 
optimal values of program (C). 

C. Simulations 

In this section, we present simulation results for the BRC 
with n = 2 destinations. We compare the multicast rate 
obtained by optimizing the relay location and the source in 
addition to relay power allocations, with the case where the 
relay is located at a naive yet seemingly interesting position: 
the centroid of triangle /\ s d 1 d 2 > an d onr y tne power allocations 











— I — Optimal relay position 














□ Centroid relay positior 
— ^— Relative gain. 











































































































4000 4500 5000 5500 

e destinations triangle (meter-square) 



Fig. 6. BRC with n = 2 destinations. Rates for r located at the centroid of 
A s d 1 d 2 > a ' th e optimal position, and relative gain. 

are optimized. The simulations are run for an increasing size 
of the area of A s d 1 d 2 an d a random network topology for each 
area is chosen. 

Figure [6] shows the maximum multicast rate (blue and red) 
for optimal and centroid relay positions respectively. The SNR 
is normalized to 1 . Note that the actual values of the rates 
are not as important because of the normalization, for higher 
power values, the rate would certainly have higher values. For 
increasing area of triangle A s d 1 d 2 > tne maximum multicast 
rate tends to drop, which is due to the constrained power and 
larger distances, but the relative gain goes up. This implies that 
for farther placed nodes the sensitivity of the relay location 
is higher and can produce significant gains in rate for the 
optimal relay location. Its clear from the results in Figure [6] 
that the centroid is not the optimal location. The rise in the 
relative gain becomes very strong due to the fact that the low- 
SNR regime is more sensitive to the location of nodes (hence, 
distances) that determine the hyperarc rates in the limit of 
disappearing SNR as opposed to e.g. in high SNR regime, 
where a displacement of ±e for the location of r would not 
effect the rate significantly. 

V. ARBITRARY CASE 
In this section, we answer the same set of questions but 

for arbitrarily placed source and destination nodes. Almost all 

concepts can be carried over to, straightforwardly. 

The steps of the pre-processing stage can be summarized as 

Input: {s, T} set of nodes with their cartesian coordinates. 

1) Compute convex hull C. 

2) Compute T~L S , switch functions (using Lemma 2). 

3) Compute the disjoint convex regions by superimposing 
all fc-order Voronoi diagrams of C. 

4) Compute H r (using relay hyperarc algorithm). 
Output: H = n s UH s - 

Once we have TL S and T~L S , we can compute O = [1, \% s \ x 
\H r |]. Ultimately, the optimization program could be stated as, 
Program (D): maximize (go(x)) 

subject to: gi(x) < 1, i g [l,k], (34) 
9j(x)<l, je[k + l,K], (35) 

where, the constraints d34i l are the posynomials constraints that 
can be transformed to convex convex constraints using GP 
and constraints d35l l are the non-posynomial constraints that 
are approximated using p-norm approximation. The objective 
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function go(x) represents the multicast rate. In program (D), 
a: is a vector of variables. 

Program (D), is an abstract representation of the actual 
program. Since, the program (D) is simply program (B) (but 
for arbitrary placement of destination nodes), the structure 
of (D) is the same as (B). The main difference is the pre- 
processing stage for the two cases, in this case which involves 
computation of k nearest neighbor nodes and superimposing 
them to form disjoint regions in C for distinct n-nearest 
neighbors. We would like to note, that this computation, 
although polynomially bounded, can be heavy. There are many 
polynomial time algorithms in the literature of computational 
geometry that solve this problem efficiently, l23ll . 

VI. RESULTS AND CONCLUSION 

A comprehensive and efficient solution is developed to 
model and answer the problem of optimal relay positioning 
so as to maximize the multicast rate from the source s to 
the destination set T in a low-SNR network. The proposed 
solution is a non-convex network optimization problem in its 
basic form that is difficult to solve. Using GP, switch functions 
and p-norm surrogate approximation we transform the problem 
to a convex approximation that can be solved using standard 
convex optimization algorithms. 

To abridge, the important contributions of this work could 
be summed up in the following words: we present a com- 
prehensive approach to determine the optimal relay position 
under the pretext of network optimization problem. Network 
topologies consisting single source, multiple destinations with 
the only intermediate node as relay are considered in complete 
generality on a 2-D plane. Using superposition coding and fre- 
quency division we construct a wireline like hypergraph. The 
low-SNR hyperarc model using superposition coding provides 
an interference free network model that is easily scalable to 
complex network topologies. Using the tools of computational 
geometry and network optimization, we presented a network 
optimization framework based solution that is intuitive and 
easy to understand. Also, we show that positioning the relay 
optimally significantly affect the network performance. 

In addition, the main causes for complexity in our approach 
are the pre-processing stage and the non-convexity arising 
from non-posynomial constraints upon GP transformation. The 
former reason could be somewhat tolerated, as generally for 
solving network planning problems heavy pre-processing is 
required. In contrast, in our case the pre-processing stage 
consists of polynomial time operations at the cost of only slight 
sub-optimality in approximation. 

The questions our work answers are just a fraction of the 
interesting questions that it opens up. An interesting direction 
would be to extend this model to general multicommodity flow 
optimization problems involving more number of relay nodes. 
On the other hand, from the computational point of view, 
an interesting question is how we can efficiently build the 
exact number of paths essentially bringing down the size of 
the network optimization program. Finding other techniques 
to model this problem could be interesting, e.g. utilizing 



geometric properties of the problem and ways to bring down 
the computational complexity. 
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